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Abstract. We present an in-depth weak lensing analysis of the 
cluster MS 1 008 .1 — 1 224 based on deep multicolor imaging ob- 
tained during the Science Verification of FORS 1 at the VLT. 
The image quality (half arcsec seeing) and depth of the VLT 
images allow the shear signal to be mapped with high signal- 
to-noise and to be traced out to 1.2 Mpc, near the edge of 
the 6'.8 x 6'.8 field of view. Using BVRI color information, as 
well as 8 1 redshifts in the field from the CNOC survey, back- 
ground galaxies can be effectively separated from cluster and 
foreground objects. PSF distorsions are found to be moderate 
across the FORS images and thus easily removed. Due to the 
small statistical errors in the mass reconstruction, this dataset 
provides a testing ground where several systematic effects (e.g. 
mass-sheet degeneracy, redshift distribution of the background 
sources, cluster galaxy contamination), which are involved in 
the weak lensing analysis, can be quantified. Several methods 
are used to remove the mass-sheet degeneracy which is found 
to dominate the systematic error budget. We measure a lower 
limit to the mass of 2.6 x 10 14 M within 1 Mpc and a 
"total" mass of 5.3 x 10 14 /ijjjj M Q by fitting a softened isother- 
mal sphere. We find the mass distribution fairly uniform, with 
no significant substructures, in agreement with the virial anal- 
ysis. The availability of the CNOC redshift data and X-ray ob- 
servations on this cluster allow a comparison of different de- 
terminations of the mass radial profile. We find the lensing and 
X-ray measurements in excellent agreement, while the mass 
derived from the virial analysis is marginally (l-2tr) in agree- 
ment at radii where both methods are reliable. This analysis 
underscores the importance of systematics in the mass determi- 
nation of clusters, particularly when such a high quality dataset 
is not available or in similar studies at higher redshifts. 
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1. Introduction 

The study of the weak lensing distortion of background galax- 
ies is a powerful tool for measuring the mass distribution of 
galaxy clusters. It has long been recognized (Tyson et al. 1984) 
that the tidal field of a cluster acts by slightly modifying the 
images of distant, background galaxies. By measuring these 
small distortions the projected mass distribution of the lensing 
cluster can be derived in a model independent fashion (Kaiser 
& Squires 1993), in contrast with other mass estimators (X- 
ray measurements of intra-cluster gas density and temperature, 
galaxy dynamics). 

To date, there are approximately 20 clusters which have 
been well studied with lensing techniques (see, e.g., Mellier 
1999 for a review). In several cases the resulting masses are 
claimed to be larger (by a factor of two) than those obtained 
using different methods. At present, it is not known if this dis- 
crepancy is due to biases and inaccuracies in the X-ray and 
dynamical estimates, or to the lensing data reduction method. 
As a matter of fact, it is not even clear whether there is a real 
discrepancy (see Wu et al. 1998; see also Allen 1998). For ex- 
ample, Boehringer et al. (1998) have shown that different mass 
estimators are in excellent agreement when very good quality 
X-ray and optical data are used. 

In this paper, we study the mass distribution of the cluster 
MS1008. 1-1224 (in the following, for simplicity, MS1008) 
using high quality imaging data obtained during the Science 
Verification of FORS 1 at the VLT. The depth and seeing qual- 
ity of these observations in four different bands (BVRI), make 
this dataset ideal for weak lensing mass reconstruction. The 
large surface density of background galaxies in these images, 
as well as the sharpness of the PSF and its small variations 
across the field, significantly improve the accuracy in measur- 
ing the distortion field from object ellipticities. As a result, high 
signal-to-noise shear maps can be obtained. This allows us to 
better investigate the effect of systematics involved in the mass 
inversion procedure, which in this case exceed statistical errors. 
Moreover, the availability of radial velocities for 65 galaxies in 
MS 1008 from the CNOC survey (Carlberg et al. 1996), as well 
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as X-ray spectroscopic observations, provide independent esti- 
mates of the total mass. 

The paper is organized as follows. In Sect. 2 we introduce 
the basic lensing relations and the notation used throughout the 
paper. The observations and the imaging dataset are briefly de- 
scribed in Sect. 3. Section 4 is devoted to the weak lensing 
analysis of the cluster, where we describe in detail the mass re- 
construction method and the results obtained. A virial analysis 
based on the CNOC redshift survey of MS 1008 is provided in 
Sect. 5. A comparison between the light and the mass distri- 
bution of MS 1008, and in particular its mass to light ratio, are 
discussed in Sect. 6. Finally, we draw the conclusions of this 
study in Sect. 7. 

An independent lensing study using the same dataset has 
recently been presented by Athreya et al. (2000); moreover, the 
depletion effect on this cluster has been recently investigated 
by Mayen & Soucail (2000, preprint). We briefly compare our 
results with these studies in the conclusions. 

Throughout the paper we adopt fl = 0.3, ft\ = 0.7, Hq = 
50 h 50 km s _1 Mpc -1 . 

2. Basic relations 

In this section we briefly review the basic lensing relations. We 
mostly follow the notation of Seitz & Schneider (1997) (see 
also Schneider et al. 1992). 

Define z& to be the redshift of the lens (in this case z& — 
0.30) and £(0) the two-dimensional projected mass distribu- 
tion of the lens, where is a two-dimensional vector represent- 
ing a direction on the sky. For a source at redshift z, we define 
the critical density, Y> c (z), which roughly describes the mini- 
mum density that a lens must have in order to produce multiple 
images from the source at redshift z. This quantity is given by 



S c (z) 



c 2 D(z) 



AwGD(z d )D(z d , z) 



for z < Zd 
otherwise . 



(1) 



Here D(z) and D(z d ,z) are the angular diameter-distance 
from the observer to an object at redshift z, and from the lens 
to the same object respectively. 

The ratio between the lens projected mass density, £(0), 
and the critical density, 'E c (z), is defined as the dimensionless 
mass density k(9,z) = £(0)/£ c (z). This quantity character- 
izes the strength of the lens for sources at redshift z; in other 
words, all the lensing observables for galaxies at redshift z de- 
pend only on k(0, z). 

In general, the source galaxies do not lie at a single red- 
shift, z, but are distributed with a some distribution p(z) (here 
we define p(z) to be normalized so that J p(z) dz = 1). This 
modifies the mean dimensionless mass density, , as 



Note that the last equality holds because T, c (z) is defined to 
be infinity for foreground galaxies. We define similarly a mean 
shear (7(0)) (see, e.g., Seitz & Schneider 1997 for a def- 
inition of the shear 7(0, z)), which also depends on the last 
integral of Eq. (||). 

In the weak lensing limit, i.e. in the limit k(0, z) <C 1 for 
all 9 and z, all of the statistical lensing observables depend on 
< K (0)) 2 and< 7 (0)) 2 . 

In this analysis, the method employed to analyze galaxy im- 
ages and infer the cluster is based on the weak lensing assump- 
tion. As shown by Kaiser (1999), such an assumption is justi- 
fied for details on the mass distribution smaller than 1/ 7 ||0||. 
For the cluster considered in this paper this limit might not be 
valid in the central region, where the presence of a small arc 
suggests that the cluster is super-critical or at least nearly crit- 
ical, and our results in this region should be interpreted with 
care. 

For simplicity, we will refer in the following to k(0) and 
7(0) as averaged quantities, (k(0)} and ( 7 (0)) . 

3. Observations 

MS 1008 is a rich galaxy cluster drawn from the Einstein 
Medium Sensitivity Survey sample (EMSS, Gioia & Lup- 
pino, 1994) at z = 0.302, with X-ray luminosity Lx[0.3 - 
3.5 keV] = 4.5 x 10 44 h^ 2 erg s -1 . MS1008 was also part of 
the CNOC Survey (Carlberg et al. 1996), whose data we will 
use extensively below. This cluster was selected for the Sci- 
ence Verification (SV) of the first instruments on the VLT-Antu 
(UT1) to provide a deep multicolor imaging in B, V, R, I, J and 
K, the FORS-ISAAC Cluster Deep Field (see Renzini, 1999). 
A detailed investigation of the mass distribution using gravita- 
tional lensing techniques was one of the main scientific goals of 
this SV Programme. Observations were carried out in January 
1999. FORS (FOcal Reducer/low dispersion Spectrograph, Ap- 
penzeller et al. 1998) was used in Standard Resolution (SR) 
imaging mode, which provides a 6'.8 x 6'.8 field of view with 
0'.'2 pixels. In order to avoid excessive bleeding and scattered 
light from two 1 1 magnitude stars in the northern part of the 
field, the FORS MOS mechanism was used in occulting mode, 
with two occulting bars on each side of the field designed to 
mask out the two bright stars at each position of the dithering 
pattern (10" step). This reduced the field of view by about 12%, 
which, as discussed below, had a moderate impact on the weak 
lensing analysis. 

We summarize in Table 1 the main parameters of the opti- 
cal imaging data. Near-IR images obtained with ISAAC, which 
cover a much smaller field of view, were not used in the present 
analysis. 

The B-band image was obtained combining all the avail- 
able frames. The VRI-band images utilized in our analysis were 
obtained coadding only the best seeing images of the entire 
set released by ESO, resulting in final images with measured 
seeing of 0.5-0.6", and depth very similar to that obtained 
coadding the whole set. The flat-fielded images were regis- 
tered and combined using procedures described in Nonino et 
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al. (1999). For the coaddition, a modified version of the driz- 
zling software originally developed for HST imaging (Fruchter 
& Hook 1999) was used. The photometric zero points are ac- 
curate at 10% level, and the uncertainty on the zero-points do 
not have a significant impact on the lensing analysis. 

4. Weak lensing analysis 

For the analysis of the weak lensing signal we used the IMCAT 
package by Kaiser (Kaiser et al. 1995, hereafter KSB) with 
some refinements, as described below. 

A single master catalog was constructed by performing ob- 
ject detection on the coadded B + V + R + 1 image. Weak lens- 
ing analysis was then carried out independently on the images 
in the four filters. In particular, for each passband, we measured 
object shapes and sizes, and performed star/galaxy separation, 
PSF correction, and mass reconstruction, as described below. 



Fig. 1. The half-light radius vs. magnitude of the objects iden- 
tified on the V band. Crosses represent objects classified as un- 
saturated stars, open squares as galaxies, and dots are discarded 
objects. 

This quantity describes the object shape regardless of its size. 
Objects with large ellipticities are highly elongated along the 
direction arg(e)/2. The quadrupole moments yield also two 
other fundamental quantities: shear and the smear polarizabil- 
ities (see below). 

Finally, object detections were visually inspected and resid- 
ual spurious sources, such as noise, cosmic rays, blends, were 
removed, and bright stars were masked out. In addition, objects 
with low S/N often have tr Q < and/or dct Q < and were 
excluded. The resulting number of objects in the four bands are 
listed in Table 2 ("good" detections). 



4.1. Object detection and shape measurements 

Objects were identified using the IMCAT peak finding algo- 
rithm which utilizes a set of images convolved with top-hat 
filters of different sizes (we used filters ranging from 0.5 to 
50 pixels). For each object the finding procedure gives, among 
other parameters, the coordinates of the center, a detection sig- 
nificance parameter, and an estimate of the object's size in pix- 
els, which we denote as r g . 

The local sky level and its gradient were measured for each 
object by computing the mode of pixels values on a circular 
annulus with inner and outer radii of 16 and 32 pixels respec- 
tively. Total fluxes and half-light radii rh were measured on 
sky subtracted images using an aperture of radius 3r g , centered 
around each object. 

Shape parameters were measured using a Gaussian weight 
function with scale length proportional to r g . The use of a 
weight function is needed in order to avoid divergences of the 
measured quadrupole moments due to the noise contribution 
from distant pixels. The complex ellipticity e was thus calcu- 
lated from the quadrupole moments Q. L j using the definition (in 
the following we will use both the complex notation and vector 
notation for quantities like the ellipticity): 

Qu - Q22 + 2iQi 2 

e = ei + ie 2 = • (3) 

Wn + W22 

Table 1. Properties of optical images. "SB limit" is the 1-cr 
surface brightness limit; "seeing" is the mean FWHM obtained 
by fitting a Moffat profile to a set of bright stars. 



Band 


Exposure 


SB limit 


Seeing 




(sec.) 


(mag arcsec" 2 ) 


(arcsec) 


B 


4950 


28.25 


0.72 


V 


1350 


27.76 


0.60 


R 


1350 


27.34 


0.52 


I 


1800 


26.48 


0.50 



4.2. Stars/galaxy separation 

Star/galaxy classification was performed independently in each 
of the four images. 

Relatively bright stars can easily be identified on the half 
light radius rh vs. magnitude plot (Fig. [l]) as a narrow verti- 
cal strip, which turns over for saturated stars. The star locus 
is clearly visible against the distribution of galaxies in the rh 
vs. magnitude plot down to V = 23 (and similarly / = 22, 
R = 22.5, B = 24). Thus, more than 60 high signal-to-noise 
(S/N) stars across the field of view (see Table 2) were used to 
model the PSF pattern and to correct galaxy ellipticities. Ob- 
jects with rh < r sta rs are mainly very faint galaxies which we 
excluded in our analysis, as the errors in the determination of 
their quadrupole moments are very large. This reduces the cat- 
alog from » 6000 ("good" detections in Table 2) to w 3000 
from which background galaxies are selected (see below). 

4.3. Identification of cluster members 

The deep multicolor photometry provided by the FORS obser- 
vations, in combination with a large number of spectroscopic 
redshifts on MS 1008 available from the CNOC survey (Yee et 

Table 2. The number of detections obtained in the four bands. 
The final number of lensing galaxies used for each band is 
given by the last line ("Relevant galaxies"). 



Bands 



Objects 


B 


V 


R 


I 


Total detections 


7969 


7969 


7969 


7969 


"Good" detections 


6369 


6315 


6213 


5395 


Stars 


66 


60 


76 


89 


Galaxies 


3042 


3429 


3584 


3178 


Cluster members 


253 


267 


269 


265 


Background galaxies 


2750 


3125 


3281 


2879 


Relevant galaxies 


1709 


1844 


2004 


1715 



CNOC galaxies 81 (58 cluster galaxies) 
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Fig. 2. The inferred PSF correction for the V image. The seg- 
ments are proportional to the correction and are oriented as the 
PSF (the longest ones corresponding to 4% ellipticities). Short 
segments outside the grid represent the uncorrected ellipticities 
of stars. 
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Fig. 3. The distribution of star ellipticities on the complex plane 
for objects identified on the V image. The "plus" signs rep- 
resent the uncorrected ellipticities, while the "multiply" signs 
represent the ellipticities corrected for the PSF. Similar plots 
are obtained for the other pass-bands. Note that all the points 
are localized near the origin of the plane, i.e. the correction is 
very small for FORS images. 



Fig. 4. Color-color diagram for the cluster. The open circles 
represent galaxies assumed to belong to the cluster; the filled 
circles represent galaxies with measured redshift in the range 
0.29-0.32; the dots are the other galaxies; crosses are residual 
galaxies with R < 21.5 (see Fig. ph. 



Fig. 5. Color-magnitude diagram for the cluster. Symbols are 
as in Fig. |j. 



al. 1998), allows the cluster galaxies to be effectively separated 
from the foreground and background galaxy population. Of the 
8 1 galaxies with measured redshifts in the FORS field, there 
are 58 cluster galaxies with 0.29 < z < 0.32, which lie along a 
distinctive "finger" of galaxies in the color-color plot (Fig. Q). 
This clump in the V — I vs. B — R includes most of the reddest 
early type galaxies in the cluster (with V — I w 1.8), as well as 
bluer cluster members, which are not simply distinguishable in 
a color-magnitude plot. Fewer late type cluster galaxies extend 
into a region dominated by foreground and background galax- 
ies (filled dots with V — I < 1.4). Thus, we are confident that 
the "finger tip" in Fig. |] encompasses the majority of the clus- 
ter galaxies (254), 23% of which are spectroscopic members. 

The color-mag diagram in Fig. |] shows the same selected 
objects, mostly lying along the cluster red sequence. A fur- 
ther magnitude cut is shown, flagging all the remaining objects 
brighter than R = 21.5 as likely cluster or foreground galax- 
ies (crosses in Fig. Q and ||). Approximately 3000 objects were 
thus selected as "background galaxies". 

4.4. Determination of the Shear 

For gravitational lensing analyses, there have been several in- 
vestigations to determine the optimal method for measuring 
galaxy shapes, correcting for non-gravitational shape distor- 
tions, and calibrating the gravitational shear estimates from 
the galaxy shape data (e.g., KSB; Luppino & Kaiser 1997; 
Kuijken 1999; Kaiser 2000). For ground based data, the stan- 
dard approach adopted in most of the literature follows the 
KSB formalism to correct for PSF anisotropies and the Lup- 
pino/Kaiser (LK) algorithm for calibrating losses due to seeing 
and pixelization. We adopted this procedure, which, qualita- 
tively, determines the shear map correcting for the isotropic and 
anisotropic part of the PSF, using a set of high S/N stars in the 
field. 

Having adopted the method for determining the shear, there 
are still several technical issues in the implementation of the 
algorithm that can introduce a bias into the results. As an at- 
tempt to minimize such effects (which we describe in more de- 
tail below), we performed a double-blind lensing analysis: two 
of the authors did independent, complete shear analyses, us- 
ing custom designed software implementations, and somewhat 
different estimators to calibrate the effects of seeing and PSF 
anisotropies. The final results from both studies were in agree- 
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ment, we thus present the results from only one implementation 
below. The second approach is described in the Appendix. 

Defining e s to be the unlensed, source ellipticity, we write 
the observed ellipticity as 

e^el+P^+P^pj, (4) 

where we have used the Einstein convention on repeated in- 
dexes. This equation simply states that, to the first order in 7 
(weak lensing limit), the observed ellipticity depends linearly 
on the shear 7. The quantity P 7 is called pre-seeing shear po- 
larizability, and is given by 

/ P sh \ 

p7 _ psh _ psm / 1 kk \ (c\ 
ij ij ij \ psm / ' yJ > 

\ r kk I * 

Here P sh and P sm are, respectively, the post-seeing shear po- 
larizability and the smear polarizability; both quantity are cal- 
culated by IMCAT using the object quadrupole moments. The 
symbol ( • }* indicates that the average is taken over the stars. 
We refer to Luppino & Kaiser (1997) for further details. 

The above relation was applied as follows. We calculated 
for each detected star the quantity P|fc/-P|j?> an( ^ fitted the 
derived values with a second order polynomial. The quantity 

/ gives information on the isotropic part of the PSF, 
and, as expected, we noticed little variation of this quantity 
across the field. 

The anisotropic part of the PSF was obtained in a simple 
way by applying Eq. (Q) to our set of stars. For stars, 7 = 
(because there is no lensing effect), and e s = (images of stars 
would be perfectly round without PSF anisotropies). In conclu- 
sion, we have p — (P sm ) e. As before, in order to obtain a 
smooth map of the anisotropies, we have fitted the p with a sec- 
ond order polynomial, p{6), which is used in Eq. (Q) to perform 
the correction. Figure || shows the fitted polynomial p(8) in the 
V band; the plots of the original and corrected ellipticities of 
stars are given in Fig. [3] (the plots for the other bands are very 
similar). 

The determination of the shear map was handled with spe- 
cial care. First, galaxies characterized by a small size, a low 
value of the detection significance v, or a small P 7 have been 
excluded from the catalog. For these galaxies errors on the de- 
termination of the shape parameters are quite large; moreover, 
they are significantly affected by the isotropic part of the PSF 
and thus they do not contribute to the lensing signal but only 
introduce noise. For the remaining galaxies we have adopted a 
weighted median defined in the following way. For each galaxy, 
we have evaluated an "individual shear estimator" given by 

7 w = ((p^w):i,,(xS (n) . (6) 

Then we assign, for each point of the field of observation 
Q, a weight W M (6) = W{6,0^) to each galaxy. Finally, 
we evaluate for each point 6 of the field a weighted median of 
both 7i(0) and 72(f)- [The weighted median x of a set of N 
real values x n with weights w n is defined as the number such 
that the sum of weights on n such that values x n < x is equal 



to the sum of weights on n such that x n > x. In other words, 
we must have 

^ w n sgn(a;„ - x) = , (7) 

n 

where sgn(a;) = x/\x\ is the function "sign" defined to be zero 
for x = 0. Since in general the condition (Q) cannot be ex- 
actly satisfied, a linear interpolation is used.] Alternatively, one 
could use a simple weighted average. Simulations, however, 
have shown that the median behaves much better. This is due 
essentially to a small number of galaxies with very poor shape 
determinations. In such a situation, clearly, the average is not 
a robust shear estimator, while the median is. As an alternative 
solution, one could use a suitable supplementary weight to take 
into account the (estimated) error on the determination of the 
shape parameters of each galaxy. This solution is not easy to 
implement because it is not always obvious if the shape of a 
galaxy has been measured correctly or not. Moreover, the use 
of supplementary weights require a fine tuning of the weighting 
method and this clearly leads to some arbitrariness on the re- 
construction. We note that a straightforward application of the 
mean (standard practice in the literature) instead of the median 
introduces in the reconstructed shear map spurious "substruc- 
tures" (i.e. secondary maxima). This is merely due to small, 
faint sources for which the measurement of the ellipticity is 
unreliable. 

For the reconstructions shown below the weight W(6, 6') 
has been implemented using a Gaussian of argument \\0 — 0'\\ 
and scale length aw = 30" (6' is the galaxy position). As 
explained by Lombardi & Bertin (1998a, 1998b), the mea- 
sured shear map has noise properties that strictly depend on the 
weight function used. In particular, we obviously cannot expect 
to detect on the shear map details smaller than the characteris- 
tic scale of the spatial weight function. For the reconstructions 
shown here we used a Gaussian spatial weight function with a 
scale length aw- 

4.5. Cluster mass 

The dimensionless projected mass map, k(8), was obtained us- 
ing an optimized reconstruction algorithm (Seitz & Schneider 
1996; Lombardi & Bertin 1998b; Lombardi 2000). 

Along with k(6), we also calculated a K (6), an estimate of 
the error on the reconstructed mass distribution. As described 
in Lombardi & Bertin (1998b), this estimate takes into account 
the intrinsic statistical error of weak lensing reconstructions, 
basically proportional to ^|e s | 2 ). We stress that a K (0) is a local 
estimate of the error; correlation on the scale of the weighting 
function is expected on the reconstructed map. 

In order to obtain the cluster mass in physical units we need 
to estimate the redshift distribution p(z) of the background 
galaxies for each passband. 

A statistical estimate of p(z) in each magnitude bin can be 
obtained by resampling the catalogue of photometric redshifts 
of Fernandez-Soto et al. (1999) in the Hubble Deep Field North 
(HDF-N). The redshift distribution of the HDF-N derived from 
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Fig. 6. The redshift distribution p{z) of the background galax- 
ies in the B, V, R, and I filters derived by resampling the catalog 
of photometric redshift in the HDF-N. 

photometric methods is believed to be fairly reliable, given the 
large number of spectroscopic redshifts available in the field 
and the good photometric accuracy of the HDF-N images. The 
drawback in using the HDF-N as a reference field is of course 
that, for such a small field, p{z) might not be representative. 
However, this is currently one of the few fields from which 
p(z) can be estimated down to our magnitude limits in an em- 
pirical fashion, without resorting to galaxy evolution models 
for N(z, m). In our analysis, we assumed that the redshift dis- 
tribution of background field galaxies, as a function of magni- 
tude, matches the one in the HDF-N. We neglected the mag- 
nification bias, i.e. the brightening of the background galaxies 
due to gravitational lensing, which should not be substantial in 
MS 1008 given the lack of prominent strong lensing features. 

If iVb (i) is the observed number of background galaxies in 
our image, in the i-th magnitude bin, and N the total number 
of galaxies detected in the image, we can write 



N h (i) 
N 



(8) 



Here p rc i(z; i) is the redshift distribution in the reference field, 
in the i-th magnitude bin, which can be written as 

1 -/ 



Pr C {(z;i) 



N Ief (i) 



>.(«) 



(9) 



N Te f(i) is the total number of galaxies in the reference field in 
the t-th magnitude bin, and z^ n ' is the (photometric) redshift of 
the 7i-th galaxy in the reference field (the sum running on all 
galaxies n found in a given magnitude bin). 



20 25 
mag 

Fig. 7. Observed number counts (solid line) for the galaxies in 
the four bands. The upper and lower solid lines correspond, re- 
spectively, to the number counts of all galaxies in the FORS 
field and those identified as background galaxies. We have su- 
perimposed (dashed line) the number counts of the photometric 
redshift catalog in the HDF-N reference field (Fernandez-Soto 
et al. 1999). The shaded areas are number counts (with statisti- 
cal errors) from the literature (Arnouts et al. 1997; Williams et 
al. 1996). 



Figure ^ shows the redshift distribution for the background 
galaxies, p(z), derived with this method. A comparison of the 
observed number counts for our field with those in HDF-N is 
given in Fig. [7| The thick solid histogram represents the num- 
ber counts of background galaxies, Nt,(i). The overall num- 
ber counts (thin solid histogram) clearly shows the excess at 
bright magnitudes due to cluster galaxies. The number counts 
of HDF-N galaxies with photometric redshifts (dashed curve) 
represent N lc f(i) in the above notation. 

The knowledge of p(z; i) = p TC f(z; i) for the background 
galaxies allows us to evaluate the projected mass distribution 
with an optimized method, which has the advantage of mini- 
mizing the errors on the mass distribution by properly weight- 
ing the contribution of each magnitude bin to the density £(0). 

For each magnitude bin i, we can define an effective crit- 
ical density as T,(z c ff^) — (E c T 1 (z)) , where the average 
is evaluated using the redshift probability distribution p(z; i) 
(bins of 1 magnitude were used). By using a weight propor- 
tional to l/J^l(z e s t i) in the average of Eq. (|), we take into 
account the fact that the lensing signal measured for galaxies 
in the bin i decreases as the effective critical density increases. 
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It is well known that the so-called mass-sheet invariance 
ultimately limits the accuracy in the determination of the total 
mass of the cluster, so that only a lower limit is strictly provided 
by the weak lensing analysis. 

A popular method to remove the mass-sheet degeneracy is 
to assume that the cluster mass density vanishes at large dis- 
tance from the cluster center. In doing this, we computed the 
radial profile of the mass map (see below) and set the minimum 
of this profile to zero. In this way, we almost certainly under- 
estimate the mass density because the field of view is not large 
enough to include distant regions of the cluster. We also note 
that taking the minimum of the two-dimensional mass map is 
generally more risky. In fact, because of the noise, the mea- 
sured mass map is not bound to be positive everywhere; on the 
other hand, the expected noise in the radial profile is signifi- 
cantly smaller (see shaded area in Fig. [ll]) and can be ignored 
for the purpose of removing the sheet degeneracy. 

Another method which is often used is to assume some par- 
ticular model for the projected cluster mass distribution, and 
to remove the sheet degeneracy by requiring that the observed 
mass matches the model. In the case of a singular isothermal 
sphere the correction to perform for the measured integral ra- 
dial profile can be written explicitly. If the mass-sheet degener- 
acy in the observed radial profile A/ b s (r) has been removed by 
assuming that the mass vanishes at a given radius, r ap , i.e. by 
setting (/i(r ap )) = 0, then the corrected mass profile M corr (r) 
is given by 

M corr (r) = ^2£M obs (r) . (10) 

Note, in particular, that when r = r ap the mass is underes- 
timated by a factor 2. In our mass estimates we have used 
both techniques described above. 

4.5.1. Projected mass distribution and radial profiles 

The mass maps obtained in the four filters are shown in Fig. [| 
in physical units. The average mass map (i.e., the arithmetic 
mean of the four independent mass maps) is also shown in 
Fig. [| superimposed on the V band image of the cluster. The 
peak of the projected mass appears to be displaced ~ 30" north 
of the cD galaxy. A significant extension of the core to the north 
is also apparent. We should emphasize that, although the mass 
reconstruction is performed independently in the four images, 
the shear maps are not completely independent because we ba- 
sically are using the same set of galaxies. In other words, we 
expect that the error on the averaged mass map is smaller than 
the error in the single bands (because we are reducing errors 
on the ellipticities measurements), but by a factor less than two 
(because we use the same galaxies). 

The radial profiles shown in Fig. |l0| were obtained by aver- 
aging the mass distributions on annuli of increasing radii. The 
same center was used for all bands, corresponding to the max- 
imum of the lensing signal. 

The azimuthal average was weighted taking into account 
the local estimate of the error for the lens mass distribution, 
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Fig. 11. The radial mass profile obtained in the B filter with 
the errors band. The dashed curve is the result of a fit with a 
smoothed isothermal sphere. Very similar results are obtained 
for the other filters. 



which depends on the local number of galaxies and on the 
smoothing function used. The use of a weighted average is cru- 
cial at larger radii, where two large regions are masked out. The 
variance cr^(6) was also used to estimate the errors on the ra- 
dial profiles. These errors should be taken with caution, since 
different points of the mass profile are correlated. The differ- 
ential radial profiles (Fig. [k| left) are in excellent agreement 
with each other, and their scatter is smaller than 5%. This can 
also be taken as a rough estimate of the error in the mass den- 
sity. The cluster mass as function of radius can be read off 
the cumulative radial profile in Fig. [l(] (right panel). In or- 
der to make these cumulative mass profiles we have broken 
the mass-sheet degeneracy by assuming, for each passband, 
that the (differential) radial profile vanishes at large radii 
(more precisely, at r = 3' 40", corresponding to our field 
of view). A comparison of different determinations of M{< r) 
in the four passbands also gives a quantitative estimate of the 
systematics involved in the measurement of the total mass, in 
particular those associated with the way the mass-sheet degen- 
eracy has been removed (which is responsible for the larger 
scatter observed in the cumulative profiles with respect to the 
differential radial profiles). For example, the mass out to lh^Q 
Mpc is in the range (2.7 3.0) 10 14 M Q , i.e. - 10%. 

The radial mass profiles can also be fitted with a softened 
isothermal model (e.g Schneider et al. 1992). This model is 
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Fig. 8. The mass distribution obtained from the four bands in units of M h^o pc -2 . The solid line represent the contour of zero 
mass density. Top Left: B band; Top Right: V band; Bottom Left: R band; Bottom Right: I band. 



Fig. 9. The mass distribution, obtained by simply averaging the mass maps in figure [| superimposed on the V image of MS 1008. 
Dashed contours are labeled in units of M Q h^o pc~ 2 . Solid contours correspond to zero mass density. The two blank regions on 
the top correspond to two bright stars masked out during the observations. 



characterized by the core radius r c and the central density So; 
the projected mass distribution is then 

(l + ||0/r c || 2 ) 3/2 



Only the central 3' of each profile were fitted with a smoothed 
version of the model Eq. (pT|), with the same smoothing length 
as the one used for the shear map (see Lombardi & Bertin 
1998b for a discussion the relationship between the smooth- 
ing scale used and the resulting mass map). The results of this 
fit are shown in Fig. [II] and summarized in Table 3. Note that 
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Fig. 10. The radial profile (left) and the cumulative radial profile (right) of the mass distributions obtained in the four bands. 



the total mass within 1 Mpc obtained from the fit is roughly 
the double of the mass assuming that the radial profile vanishes 
at large radii. This discrepancy is related to the mass-sheet de- 
generacy and clearly represents a serious limitation for a "total" 
mass estimate. We stress in particular that the correction to be 
applied to the "detected" mass in order to obtain the total mass 
strongly depends on the radial profile assumed for the cluster 
(the correction is exactly 2 for a singular isothermal sphere, see 
Eq.©). 

4.5.2. Integrals of the shear 

An alternative method to evaluate the radial profile is to use in- 
tegrals of the shear (Fahlmann et al. 1994; Squires & Kaiser 
1996). In particular, for a point located at polar coordinates 
(r; ip) we define the tangential shear as 7 t — —(71 cos2</? + 
72 sin 2<p). An estimator of the mass profile is given by 

C(0i) = (k)(0 1 )-(k)(0 1 ,6 2 ) 



Table 3. Best fit parameters obtained in the B band. 



Parameter 



Value 



r c (arcsec) 
So 

ID velocity dispersion 
Mass(r < 1 Mpc) (fit) 
Mass(r < 1 Mpc) (direct) 



35 
265 
850 
5.3 x 10 14 
2.6 x 10 14 



arcsec 
M hso pc" 
km s _1 

M 
h7j M 



*0--9ilH)Jei Jo 



dr 



271 



lt(r',(p) 



dip 



(12) 



which measures the mean dimensionless mass density (k)(6i) 
inside a disk of radius 9\ relative to (k)(0i, 62), the mean in an 
annulus of radii 9% and 9 2 . 

We show C(9i) in Fig. [j~2| for each passband; we used #2 = 
3' 40" which defines the edge of the images. The associated 
integral mass estimator, i.e. the quantity 



Z(0i) = E c 7rC(0i,0 2 ; 



(13) 



is shown in Fig. [f2| In Fig. [f2] we also show in gray the same 
estimate as obtained independently by one of us (G. Squires) 
using a slightly different technique (see Appendix). 



5. Mass profile from the virial analysis 

The analysis of the internal dynamics of MS 1008, as traced by 
its member galaxies, has already been discussed in the liter- 
ature by the CNOC collaboration. The internal (line-of-sight) 
velocity dispersion for this cluster was originally measured by 
Carlberg et al. (1996), who found a v = (1054 ± 107) km s" 1 , 
based on the explicit background subtraction of the interloper 
contamination. More recently, Borgani et al. (2000) consis- 
tently found a v = (1042 ± 110)kms _1 , based on the com- 
pletely different method for removing interlopers, as described 
by Girardi et al. (1998, G98 hereafter). Carlberg et al. (1996) 
also estimated M200 — 2 • 10 15 M Q for the mass within the 
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Fig. 12. The (differential) mass profile estimator C, (left) and its integral, Z = M (< r) (right), for the four passbands. The shaded 
are shows the result from the weak lensing analysis with a different method (see text). 



radius ?'20o — 2.9 h^ Q Mpc, which is defined as the radius en- 
compassing an average overdensity of 200 (all the scales from 
the original papers have been suitably converted to those ap- 
propriate for the reference cosmological model assumed in our 
analysis). Lewis et al. (1999) estimated the mass within a ra- 
dius r ~ 0.8 Mpc from ROSAT HRI X-ray data and found 
Mx ~ 3.6 • 10 14 h^Q M©. They also found this value to be con- 
sistent with that obtained by extrapolating to this smaller radius 
the virial mass estimate at r 2 oo with the NFW-like (Navarro et 
al. 1996) profile, fitted by Carlberg et al. (1997) for the average 
density profiles of member galaxies for CNOC clusters. 

Here, we use the galaxy data (redshift and positions) as 
provided by Yee et al. (1998) to determine the mass-profile of 
MS 1008 from the virial analysis. The procedure of the selec- 
tion of member galaxies is the same of G98. In addition, we 
also reject emission line galaxies, since they are believed not to 
be fair tracers of the cluster internal dynamics (e.g., Biviano et 
al. 1997). From G98 we also adopt the procedure for estimating 
the velocity dispersion a v and the virial mass M v i r , for which 
we provide here only a short description. As for the velocity 
dispersion, we find a v = 1033^105 ^ m s_ 1 > tnus ver y c l° se to 
previous estimates. 

In the case of a spherical system and under the assump- 
tion that mass follows galaxy distribution, it is possible to show 
(e.g., Limber & Mathews 1960; The & White 1986) that 



37T (T^i?PV 

~2 G 



C . 



(14) 



Here C is a term of surface pressure correction, which takes 
into account the fact that the radial component of the velocity 
dispersion does not vanish within the observational aperture ra- 
dius (Carlberg et al. 1997). In the following, we assume a 20% 
correction at the virialization scale, which is typical for nearby 
and CNOC clusters (G98; Carlberg et al. 1996). The value of 
the projected virial radius i?py depends on the relative (pro- 
jected) positions of member galaxies within the sampled clus- 
ter region. 

Equation ( |l4| ) provides a reliable estimate of the virial mass 
only if galaxy redshifts are available over the whole virial- 
ized region of the cluster. Unfortunately this is not the case 
for MS 1008 which is sampled only within about half of the 
virialized region (cf. Carlberg et al. 1996). While this is not a 
serious problem for the estimate of a v , whose profile flattens 
well within the aperture radius (Borgani et al. 2000), it may 
represent a limitation for Rpy, whose value increases with in- 
creasing clustercentric distance. However, at each radius, an 
estimate of Rpy can be obtained by assuming a shape for the 
density profile. For instance, G98 estimated Rpy by assuming 
a King-like profile, 

[1 + (r/r c y\ 

with core radius r c ~ 0.1 Mpc and shape a ~ 1.2 [cf. their 
Eq. (13)], which is the best fit to the average galaxy density 
profile for an extended sample of local clusters. They also veri- 
fied this estimate of Rpy to be reliable for clusters sampled out 
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to the virialization scale. In this way, we find a total virial mass 
which fully agrees with that given by Carlberg et al. (1996). 
Diamonds in Fig. [l3| show the projected mass profile obtained 
from the King-like profile of Eq. (|T5|). 

Figure 13 also shows with the triangles the projected mass 
profiles obtained by scaling the original CNOC estimate with 
the NFW-like profile fitted by Carlberg et al. (1997). Since 
no definitive conclusion has been reached at present about the 
most appropriate profile describing the cluster galaxy distribu- 
tion (see, e.g., Merritt & Tremblay 1994; Adami et al. 1998), 
we prefer here to show the results for both such profiles. Al- 
thouh the two virial mass determinations coincide on scales 
comparable to the cluster virialized regions, the choice for the 
density profile affects the mass reconstruction at small radii. 

As a remarkable result, we find an almost perfect concor- 
dance between our weak lensing mass and the X-ray mass by 
Lewis et al. (1999). We stress that the weak lensing profiles 
shown in Fig. [D]have been calculated by assuming a singu- 
lar isothermal sphere for the cluster mass distribution (we 
recall that this assumption is important when removing the 
mass-sheet degeneracy). As for the comparison with the virial 
mass reconstruction, a fairly good agreement is found at least 
for 0.4/igQ 1 Mpc, the virial mass being only slightly larger 
(at about la confidence level) than that provided by the weak 
lensing reconstruction. The difference becomes significant only 
at smaller radii, where, however, both the virial and the lensing 
masses become more uncertain. 

These results indicate that the weak lensing mass mass re- 
construction do provide results which agrees with other meth- 
ods, based on the assumption of dynamical equilibrium, at least 
in cases when high-quality imaging data are available for a 
fairly relaxed cluster, like MS 1008 (see also Allen 1998). 



6. Mass to light ratio 

To evaluate the mass to light ratio, we constructed a smoothed 
map of the light distribution in the B band, using a Gaussian 
kernel with the same scale as the one used for the shear map. 
The light density distribution is shown in Fig. [l4| with overlaid 
mass contours, and the cluster galaxies used in our analysis. 
Galaxy fluxes were measured with the SExtractor package 
(Bertin & Arnouts 1996), which ensures better object deblend- 
ing in the cluster core with respect to IMC AT. 

Luminosities were computed by assuming a k-correction of 
0.7 magnitude in V band (Fukugita et al. 1995), appropriate for 
ellipticals which dominate the cluster galaxy population. Lumi- 
nosity evolution of early type galaxies can be neglected at this 
low redshift. In our cosmological model the distance modulus 
is DM = 51og[D(z d )(l + z d ) 2 /10pc] = 41.77 -5 log h 50 . 
The resulting rest-frame V-band luminosity summed over all 
the cluster galaxies is Ly = 5.6 x 10 12 /i^~ 2 Lq. This esti- 
mate includes galaxies as faint as M* — 4, hence the correc- 
tion for incompleteness should be negligible. Carlberg et al. 

8.8 x I0 12 h^ L Q (g = 0.1) in the Gunn r 
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(^50 Mpc) 



1.0 



found L 2 oo 

band, at r = i? 20 o ^ 2.3^ Mpc. Using (V - R) Ea = 0.9 
and (r — R)eii = 0.3, our estimate translates to LQ unnr ~ 
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Fig. 13. Comparison of different projected mass estimates. Dia- 
monds indicate results based on the virial mass from the present 
analysis, rescaled with the King-like density profile of Eq. (p^|). 
Triangles are for the virial mass by Carlberg et al. (1996), 
rescaled with the NFW-like profile. Errorbars are la uncertain- 
ties propagated from the statistical errors in the a v estimate. 
The solid square is obtained by projecting the X-ray mass of 
Lewis et al. (1999), under the assumptions of spherical symme- 
try and NFW-like profile. The two curves are the mass profiles 
derived from weak lensing in two different bands. 

8.4 x 10 12 ft^ 2 L at r ~ 1.4/igg 1 Mpc for q = 0.1; hence it is 
in very good agreement with L200 estimated by Carlberg et al. 

Figure [14] clearly shows that, as commonly observed in 
clusters (e.g. Clowe et al. 1998), the light roughly traces the 
mass. The NE elongated morphology, clearly visible in the 
smoothed light distribution, is also shared by the intra-cluster 
gas, as visible in Fig. [l4|, which shows the ROSAT-HRI X-ray 
contours overlaid on the mass map (grey scale). Interestingly, 
the X-ray centroid is also displaced 15" North of the cD, i.e. 
nearly coincident with the peak of the weak lensing signal. The 
X-ray HRI data are not deep enough to probe the outer regions 
of the cluster. We also find some evidence for a drop of M /L 
in the center possibly due to the presence of the over-luminous 
cD galaxy. 

In order to properly take into account the sheet degeneracy 
in the evaluation of the mass-to-light ratio, we used an estima- 
tor similar to £ for the light distribution: 



Cv(0i) = (Lv)(0i) - (LvKh 



(16) 



where (Ly}(9) is V band luminosity averaged over a disk of 
radius 9. Thus, C(^i)/Cv(#i) yields an unbiased estimate of the 
mass-to-light ratio (assuming that the light traces the mass). 
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Fig. 14. Left: Light distribution of MS1008 (gray scale) with overlaid mass contours, as in Fig.|| Dots represent galaxies selected 
as cluster members on the basis of their color (filled dots) or redshift (crossed dots). Right: The ROSAT HRI X-ray contours 
overlaid on the weak lensing mass distribution (gray scale). 

Using this method we obtained M / L\ — 95 /150 M© /L© 
within a radius of 3', which compares favorably with the 
value found by Hoekstra et al. (1998), M/L v = (90 ± 
13)ft. 5O M /L , for the CL 1358+62, with redshift (z = 
0.33) and X-ray luminosity (i x [0.3 - 3.5 keV] = 5.5 x 
10 44 h^ 2 erg s" 1 ) very similar to MS 1008. 



7. Conclusions 

We have studied the weak lensing mass distribution of the clus- 
ter MS 1008 using multicolor imaging data obtained during the 
Science Verification of FORS 1 on the VLT. The depth, angu- 
lar resolution, and image quality across the entire field of view 
make these images a rare and ideal dataset for weak lensing 
analysis. In addition, the information on radial velocities of 
a large set of cluster galaxies provided by the CNOC survey 
allows the mass of this cluster to be derived with dynamical 
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methods and compared with that reconstructed from the weak 
lensing. These combined datasets allowed us to conduct a de- 
tailed study of the systematic effects involved in cluster mass 
estimates. 

We have found PSF distortions of FORS1 to be moder- 
ate and hence easily removed with a low order polynomial fit. 
The B through I multicolor information, as well as the red- 
shift information from the CNOC catalog, have allowed us 
to efficiently separate the background galaxies from the clus- 
ter and foreground populations. Approximately 8000 objects 
have been detected, corresponding to 200 sources arcmin~ 2 . 
Further selection lead to approximatively 1700 background 
galaxies with high signal-to-noise ratio which have been used 
for the determination of lensing shear field. We have de- 
tected a weak lensing signal out to r ~ 1.2h§ Q Mpc. The 
projected mass within l/ijj" Mpc is measured in the range 
(2.5 -7- 3.2) lO 14 ^ 1 M in the four filters, with typical statis- 
tical errors (deriving from uncertainties in galaxy ellipticities) 
of 5%. 

We have discussed the impact of systematics in the weak 
lensing reconstruction of the mass distribution in physical 
units, which can be summarized as follows. The removal of 
the mass-sheet degeneracy is inevitably model-dependent, in- 
ducing mass variations of ~ 15% at r w I^sq 1 Mpc. [Note that 
the estimate of this error is difficult and somewhat arbitrary 
since the error clearly depend on the "range" of possible mod- 
els allowed. As a result, the estimate of the total error, given 
below, is also bound to be inaccurate. We stress, however, that 
a similar problem exists for X-ray or virial mass estimate (cf., 
e.g., the different predictions using King-like or NFW profiles 
if Fig. |l3[).] By polluting the background galaxy population 
with a significant fraction of cluster galaxies, the total mass 
is biased low by ^ 5%. However, different selections of back- 
ground galaxies, as well as edge effects due to masked portions 
of the images, do affect the morphology of the shear maps. We 
have not found any evidence for a substructure of the cluster 
core; this in part is due to lack of resolution in the mass re- 
construction, set by the smoothing scale used. However, we 
have found the central region of the cluster to be elongated NS, 
with some evidence of a displacement of the projected mass 
peak respect to the cD galaxy, similarly to the X-ray emis- 
sion. Cosmic variance is expected to affect the assumed red- 
shift distribution of the background field galaxies, and hence 
the derived cluster mass. From a comparison between the Hub- 
ble Deep Field North and South, cosmic variance would intro- 
duce a systematic error of approximatively 10% in the mass. 
This can also be taken as a rough estimate of the error due 
to cosmic variance. Note also the main contribution to this 
error is from large scale structures rather than from Pois- 
son noise in the HDF-N and in our field. In other words, a 
simple error estimate which assumes that galaxies are un- 
corrected would lead to an underestimate of the error. Al- 
together, the cluster light distribution traces the mass distribu- 
tion remarkably well. We have measured a mass-to-light ratio 
of M/L v = (95 ± 22) h 50 M Q /L Q at radii > 0.5h^ Mpc. 



We have also compared the weak lensing mass profile 
with that derived from a virial analysis of the CNOC red- 
shift data, as well as with the X-ray mass. The mass de- 
rived from X-ray observations (Lewis et al. 1999) is found in 
excellent agreement with our weak lensing determination at 
r ~ 0.8/igQ 1 Mpc. Different approaches to estimate the cluster 
mass at the virialization scale produce very similar results. At 
scales 0.1 ^ h^r ^ 1 Mpc, lensing and virial mass determina- 
tions can be compared. The virial mass reconstruction at these 
small radii critically depends on the assumed density profile. 
Using a King-like profile describing the galaxy distribution of 
local clusters (G98; cf. also Eq. (|l5|)), we have obtained a mass 
which is more than 2er away from the X-ray and weak-lensing 
mass. Assuming instead a NFW profile, as fitted by Carlberg 
et al. (1997) to CNOC clusters, a much better agreement is ob- 
tained over the entire overlapping region. Note that the discrep- 
ancy at small radii (r < 1') between lensing and virial estimate 
can be explained by recalling that the lensing determination is 
affected by a Gaussian smoothing of 30", as well as a departure 
from the weak lensing approximation. The former effect leads 
to underestimate the mass in the core, while the latter will gen- 
erally bias the central mass high. Both effects clearly vanish at 
large radii. In addition, we notice that the weak lensing data 
alone cannot discriminate between different mass profiles (e.g. 
NFW versus isothermal model). 

This analysis shows that the combination of depth and good 
imaging quality provided by VLT/FORS allows us to obtain 
high S/N mass maps via optimized weak lensing reconstruc- 
tion methods. On the other hand, these data have the virtue of 
revealing important systematic errors in weak lensing analyses. 
We find that the mass-sheet degeneracy dominates the budget 
of systematic errors and ultimately makes mass measurements 
via weak lensing model dependent. 

The results obtained here are, at a first look, in good agree- 
ment with the mass reconstruction performed by Athreya et 
al. (2000). However, a detailed analysis of their study reveals 
that a different recipe to remove the mass-sheet degeneracy was 
used. Specifically, they set to zero the mass at large radii (3'.2- 
3'.6). Had they used our (model-dependent) method to break the 
mass-sheet degeneracy, their estimate would have been a factor 
of two higher. Therefore, there is a discrepancy of the same fac- 
tor between the two analysis which, at present, we are unable 
to explain. 



8. Acknowledgments 

The authors would like to thank the UT1 Science Verification 
team at ESO for the preparation of this unique data set and 
its release to the community. We are particularly grateful to 
the VLT/UT1 Commissioning Team for carrying out the ob- 
servations. We also thank Tom Broadhurst, Tomas Erben, and 
Peter Schneider for useful discussions. Support for this work 
was provided by NASA through Hubble Fellowship Grant No. 
HF-01 1 14.01-98A from the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities for 



14 



Lombardi M. et al.: Weak lensing mass reconstruction of MS 1008.1 — 1224 



Research in Astronomy, Incorporated, under NASA Contract 
NAS5-26555. 



Appendix A: Alternative Shear Analysis 

In the second analysis of the data, we followed the procedure 
used by Hoekstra et al. (1998). The anisotropic part of the PSF 
is removed on an object-by-object basis by correcting each ob- 
ject's ellipticity via 



psm 



PI 



(A.l) 



where starred quantities refer to parameters measured for stel- 
lar images. In this analysis, we calculated the stellar ellipticities 
and polarizability employing a matched radius for the Gaussian 
weight function as for the object in question. We fit a second 
order polynomial to account for spatial variations in the PSF. 

We next determine the shear from the anisotropy-corrected 
galaxy polarizations. Both the (now effectively isotropic) PSF, 
and the circular weight function, tend to make objects rounder. 
These effects may be corrected for using the 'pre-seeing shear 
polarizability' 



pi _ psh, _ psh. J psm x p 



(A.2) 



introduced by Luppino & Kaiser (1997), where the * subscript 
again refers to quantities determined from the objects classified 
as stars. 

Operationally, to determine (P 1 ), we split the faint galaxy 
catalogs into six bins based on size r g , with spacing dr g = 0.25 
pixels. We reject galaxies with very large, small, or negative 
values of P sh = (Pff + P$)/2, and P srn = (Pff + P 2 T)/ 2 
(the off-diagonal terms typically being very small). The ratio- 
nale for this cut is that the calculation of the shear and smear 
polarizabilities for these objects failed, typically either because 
the objects themselves were very faint and near the detection 
threshold, or due to the presence of very nearby object con- 
taminating the aperture over which the calculations were per- 
formed. 

In each r g bin, we calculate the median galaxy polarizabil- 
ities, P sh and p sm , and reanalyze the stellar catalog, resetting 
the value of the stellar r g to the bin size, and analyzing with a 
Gaussian with scale r g and an aperture three times this radius. 
We finally estimate the shear for each object using 



7a 



<P7) 



(A.3) 
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